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ABSTRACT 

In our previous paper, we have reported the detection of a Balmer edge absorp- 
tion feature in the polarized flux of one quasar (Ton 202). We have now found similar 
Balmer edge features in the polarized flux of four more quasars (4C09.72, 3C95, B2 
1208+32, 3C323.1), and possibly a few more, out of 14 newly observed with the VLT 
and Keck telescopes. In addition, we also re-observed Ton 202, but we did not detect 
such a dramatic feature, apparently due to polarization variability (the two observa- 
tions are one-year apart). The polarization measurements of some quasars are affected 
by an interstellar polarization in our Galaxy, but the measurements have been cor- 
rected for this effect reasonably well. 

Since the broad emission lines are essentially unpolarized and the polarization is 
confined only to the continuum in the five quasars including Ton 202 in both epochs, 
the polarized flux is considered to originate interior to the broad emission line region. 
The Balmer edge feature seen in the polarized flux is most simply interpreted as an 
intrinsic spectral feature of the quasar UV/optical continuum, or the "Big Blue Bump" 
emission. In this case, the edge feature seen in absorption indeed indicates the thermal 
and optically-thick nature of the continuum emitted. However, we also discuss other 
possible interpretations. 

Key words: quasars - galaxies: active - accretion - polarization - radiation mecha- 
nisms: general 



1 INTRODUCTION 

We have recently reported our first detection of a Balmer 
edge absorption featu re in the polarized flux spectrum of 
one quasar (Ton 202; iKishimoto. Antonucci fc BlaedEooll 
hereafter Paper I). In the present paper, we report that a 
similar Balmer edge feature is seen in the relatively high S/N 
polarized flux spectrum of four more quasars, taken with the 
Very Large Telescope (VLT) and Keck Telescope. We also 
present new spectropolarimetric data for 10 more quasars of 
a mostly moderate S/N, also taken with the VLT and Keck. 

These observations are a part of our on-going effort to 
find spectral features intrinsic to the UV/optical continuum 
of quasars, which is often called the Big Blue Bump (BBB). 
This emission component is crucial in the sense that it dom- 
inates the radiative out put of quasars, but its nature has not 
been well understood llKoratkar fc Blaesl Il999l : lAntonuccl 
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Il988tll999l) . One of the most critical problems has been the 
observed apparent lack of continuum edge spectral features 
from the emission source, which is often assumed to be an 
accretion disk. 

At least at the Lyman edge, this problem may be miti- 
gated by smearing from relativistic Doppler shifts and grav- 
itational redshifts, which can severely blur the spectral fea- 
tures produced in the local atmospheres at different radii of 
the disk (e.g. Hubeny et al. 2000). While the coolest disk 
models still exhibit a broadened absorption feature, hot- 
ter disk models generally exhibit a change in slope in the 
vicinity of the Lyman limit. This arises from a combina- 
tion of smeared out absorption and emission edges on top of 
a general rollover of the underlying continuum. Such slope 
changes might be consistent with the spectral breaks that 
are observed in the vicinity of the Lyman limit in some ob- 
jects such as 3C 273 (Kriss et al. 1999) or in composite 
spectra (Zheng et al. 1997; Telfer et al. 2002), but quanti- 
tative fits to such spectra have so far been unsatisfactory. 
In particular, detailed fits to the spectral break observed in 
3C 273 still produce a local emission bump that is not in 
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agreement with the data (Blaes et al. 2001). This bump in 
the best fit disk model arises from Lyman emission edges 
formed in the atmospheres between ten and thirty gravita- 
tional radii. Relativistic smearing of this bump is minimized 
because 3C 273 is a superluminal source and the putative 
disk is therefore constrained to be nearly face-on to the line 
of sight. More generally, disk models that provide qualitative 
fits to the Lyman spectral break still do not provide enough 
emission in the far ultraviolet and soft X-rays to fit the com- 
posite spectra because the thermal continuum rolls over too 
soon (Laor et al. 1997; Blaes 2004). Some additional source 
of emission, e.g. Comptonization (Czerny & Elvis 1987; Lee, 
Kriss & Davidsen 1992), is probably required. 

Compared to the Lyman edge region, accretion disk 
models place the main emitting source for the Balmer edge 
region considerably farther out in the gravitational poten- 
tial well. In detail, such models predict substantial contri- 
butions to both the Lyman and Balmer edge regions from 
broad ranges of overlapping radii, but the Balmer edge still 
has contributions from considerably larger radii, and there- 
fore will be subject to less smearing. (There are actually 
considerable technical problems in making concrete spectral 
predictions of the Balmer edge from accretion disk mod- 
els because of the prese nce of hydrogen ioni zation zones at 
these larger radii - see Hu benv et al] 120001 for discussion. 
It is noteworthy, however, that the Balmer edges in these 
models is always in absorption.) In addition, unlike the Ly- 
man edge, the Balmer edge is not a resonant feature, so the 
confusion with foreground absorption is much less of a prob- 
lem. Thus, although detailed spectral predictions from disk 
atmosphere models of the Balmer edge are yet to come, tar- 
geting the Balmer edge may well be advantageous compared 
to the Lyman edge. 

However, the observation of the Balmer edge intrinsic 
to the BBB emission has been impossible, since it is es- 
sentially buried under the Balmer continuum, high-order 
Balmer emission lines, and Fell emission lines (collectively 
referred to as the small blue bump) which are thought to 
be coming from the broad emission line region (BLR). The 
main immediate goal of our effort is to overcome this ob- 
stacle by looking at the polarized flux of the quasars with 
a small polarization (P ~ 1%). Note that these are normal 
quasars in terms of their spectral energy distribution, and 
not classified as blazars. About 1/3 of the normal quasars 
are k nown to be polarized with P > 0.6% l|Berriman et alJ 
1990). Among them, at least in some cases (and possibly in 
many cases), the polarization is found to be confined only to 
the continuu m, and the broa d emission lines are essentially 
unpolariz ed <IAntonuccil988[. show ing the data of Miller and 
Goodrich: ISchTnTh^^^mit,hl200nl) . In these cases, the polar- 
ized flux is considered to arise interior to the BLR, and thus 
it is likely to show the intrinsic Balmer edge feature of the 
BBB emission with all contaminating emission from outside 
the nucleus scraped off and removed. 

Detailed investigation of the polarized flux requires a 
long integration with a large telescope. In our first such deep 
integration on one quasar with the Keck telescope, we have 
indeed found a Balmer edge feature in the polarized flux of 
the quasar. Here, we report our spectropolarimetric observa- 
tions of 14 more quasars with the VLT and Keck, presenting 
a similar Balmer edge feature in four more quasars. We also 



Table 1. List of symbols and acronyms used in this paper. 



symbol/ definition 
acronym 

I Stokes / parameter, or the total flux F. 

Q, U Unnormalized Stokes parameters in the sky co- 

ordinates, i.e. having the reference axis along 
the North-South direction. 

q,u Normalized Stokes parameters, q = Q/ 1 , u = 

u/i. 

Q' ,U' ,q' ,u' Rotated Stokes parameters with the reference 
axis taken to be along the PA of the polariza- 
tion of the object. Thus, Q' (or q' X F) is essen- 
tially equivalent to the polarized flux P X F. 

rp Ratio of polarization at 2891-3600A to that at 

4000-4731A. 

APA Difference between the PA of the optical po- 

larization and that of the radio structural axis. 

A max The wavelength at which an interstellar polar- 

ization is maximum for a given line of sight. 

P max The magnitude of an interstellar polarization 

at Amax • 

PA Position Angle. 

BLR Broad Line Region. 

BBB Big Blue Bump. 

ISP Interstellar Polarization. 



discuss the re-observation of Ton202. Some of the prelimi- 
nary results have been published in Kishimoto et al. (2004) . 

Finally, note that the polarization P and the polarized 
flux P x F are different quantities which must not be con- 
fused. It has been previously established that P tends to be 
low in the Balmer continuum region, but our work, which 
has important new physical implications, is that PxFin the 
same wavelength region is also well below the extrapolation 
from longer wavelengths in many cases. 

For convenience, we summarize the symbols and 
acronyms used throughout this paper in Tabled 

2 OBSERVATION 
2.1 VLT observations 

We observed 11 quasars on 13 and 15 September 2002 
(UT) using Focal Reducer/Low Dispersion Spectrograph 
(FORS1) on the Unit Telescope 3 (Melipal) of the Very 
Large Telescope (VLT). These observations are summarized 
in Table [5] The observation of each target were done in 
one or more sequences with one sequence consisting of four 
frames with four different waveplate positions. The expo- 
sure time in Table [5] refers to the exposure time for one 
sequence multiplied by the number of repetitions of the 
same sequence. The grism 300V was used at a dispersion 
of 2.6A per pixel. The seeing varied between ~ 0."5 and 
~ 1" over the nights: the slit width was set to l."5 giving 
a spectral resolution of 20A, except for some frames taken 
in the first half of 13 Sep (see Table |5J when the seeing was 
mostly 0."5 — 0."7 and the slit width was set to 1" (res- 
olution 13A). The spatial sampling was 0."200 per pixel. 
The slit direction was fixed north-south, and the instrument 
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Table 2. Log of VLT/FORS1 observations in 2002. The slit PA was fixed at 0°. 



name 


other 


z 


exp. 


date 


slit 


filter 




name 




(min 


) 






(") 




CinnnQ_i_i f^c 
v^UUUo-t-loo 




U.40U 


20 


X 


1 


13 


Sep 


1.0 


OLtO / 








20 


X 


1 


15 


Sep 


1.0 


OLtO / 






u. roo 


10 


X 


1 


13 


Sep 


1 ,o 


none 






U.010 


.30 


X 


3 


13 


Sep 


1 ,o 










16 


X 


1 


13 


Sep 


1 .O 










30 


X 


2 


15 


Sep 


1 ,0 










30 


X 




15 


Sep 


1 ,o 


in 'lie 








20 


X 


1 


15 


Sep 


1.5 


none 


Q0350-073 


3C94 


0.962 


20 


X 


1 


13 


Sep 


1.5 


GG435 


Q0405-f23 


PKS 


0.573 


10 


X 


1 


13 


Sep 


1.0 


GG435 


Q04f4-060 


3C110 


0.775 


10 


X 


1 


13 


Sep 


1.5 


GG435 


Qf9f2-550 


PKS 


0.402 


20 


X 


1 


13 


Sep 


1.0 


GG375 








20 


X 


5 


15 


Sep 


1.5 


GG375 








30 


X 


1 


15 


Sep 


1.5 


GG375 


Q2fl5-305 


PKS 


0.979 


10 


X 


1 


13 


Sep 


1.0 


GG435 








10 


X 


3 


15 


Sep 


1.0 


GG435 


Q225f+lf3 


4C11.72 


0.326 


10 


X 


1 


13 


Sep 


1.0 


none 


Q2308+098 


4C09.72 


0.433 


30 


X 


2 


13 


Sep 


1.0 


GG375 








30 


X 


2 


13 


Sep 


1.5 


GG375 


Q2349-014 


4C-0f.6f 


0.474 


10 


X 


1 


13 


Sep 


1.0 


none 



Table 3. Log of Keck observations on 4 May 2003. The slit width was l."5. 



name 


other 


z 


exp. 


slit PA 


grism 




name 




(min) 


(°) 




Q1004+130 


4C13.41 


0.240 


30x1 


58 


300/5000 


Q1208+322 


B2 1208+32 


0.388 


30x3 


30 


300/5000 








30x1 


30 


400/3400 


Q1425+267 


Ton 202 


0.366 


30x2 


143 


300/5000 








30x1 


97 


400/3400 


Q1545+210 


3C323.1 


0.264 


30x2 


85 


300/5000 








16x1 


85 


400/3400 



was equipped with an atmospheric dispersion compensator 
iAvila. Rupprecht. fc Beckerslll997l) . 

The CCD frames were bias-subtracted using an aver- 
aged bias frame, and flat-fielded using internal screen flats. 
Obvious cosmic-ray hits were removed and interpolated us- 
ing neighboring pixels. The o-ray and e-ray spectra were 
extracted with 2. "2 windows. Then these were combined 
to produce normalized Stokes parameters q and u, and the 
Stokes parameter /, following iMiller. Robinson fc Goodrich! 
(1988). The zero point of the position angle (PA) of the 
polarization was determined using the polarized standard 
stars Hiltner 652, BD+25 727, and BD-12-5133, with an es- 
timated accuracy of 0.3°. The instrumental polarization was 
checked with the unpolarized standard stars HD 13588 and 
HD42078 and estimated to be less than ~ 0.05%. The flux 
was calibrated using the flux standard star CD 50. 

In many of our observations, we have used two different 
order-sorting filters to avoid second-order contamination; we 
have chosen for each object the filter which blocks the sec- 
ond order spectrum at the H/3 line position in the observed 
frame. We observed one of the objects (3C95) also without 
the filter to cover the near-UV wavelengths. 



2.2 Keck observations 

Four quasars listed in Table El were observed on 4 May 2003 
(UT) with the L ow Resolution Imaging Spectrograph (LRIS; 
lOke et aljfl995h on the Keck-I telescope, with a new blue- 
sensitive CCD (commissioned on 4 June 2002, after the ob- 
servations in Paper I) on the blue arm of the spectrograph. 
Two grisms, 300/5000 and 400/3400, were used with dis- 
persions of 1.4A per pixel and l.oA per pixel, respectively. 
The slit width was l."5, giving a spectral resolution of ~14A 
with grism 300/5000 and ~10A with 400/3400. The spatial 
sampling along the slit was 0."135 per pixel, and the data 
were obtained with the slit roughly at the parallactic angle. 
The seeing size was estimated to be about 0."8 in the first 
half of the night. We probably were affected by some cirrus 
especially in the second half of the night. 

The 2D spectral images were reduced in the standard 
manner. The bias level was subtracted using bias frames and 
overscan regions. The o-ray and e-ray spectra were extracted 
with 2. "8 window. Then these spectra t aken at 4 posi t ions o f 
the waveplate were combined following IMiller et al.l lll988}) . 
The zeropoint of the PA was calibrated with the polarized 
standard star HD183143. The instrumental polarization was 
checked with the unpolarized standard stars HD94851 and 
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BD+28 4211 and also with the observations of Ton202 with 
two different slit PAs. The former two observations indicated 
possible instrumental polarization of P ~ 0.13% roughly in- 
dependent of wavelengths at a certain PA in the instrumen- 
tal coordinates. However, the latter observations of Ton202 
at two PAs showed no systematic difference (the correction 
with 0.13% results in systematic difference of this level). 
Therefore, we did not correct for any instrumental polariza- 
tion. The flux was calibrated using the standard stars Feige 
34 and BD+28 4211. 

For the observations with the grism 300/5000, there was 
no order sorting filter available, and the new blue arm CCD 
is very blue sensitive. Therefore, the spectra at the red side 
are contaminated by second-order light. To correct for the 
contamination in flux, we implemented short exposures on 
Ton202 and B2 1208+32 with dichroic D460, which cuts out 
the light with wavelengths longer than ~ 4600A. This pro- 
vided the count rate spectrum for the second order light (the 
spectrum < 5500A was masked). We scaled and subtracted 
it from the Stokes I spectra of these two objects. Therefore, 
the total flux spectrum is expected to contain very little 
contamination. However, the polarization measurement was 
not corrected (the observations with D460 were done only 
at one position angle of the waveplate), so the polarized 
flux spectrum could be slightly affected. For 3C323.1 and 
Q1004+130, the spectra were not corrected, but these have 
lower redshift than the other two objects, so the wavelength 
range up to H/3 line is probably not affected severely. 



2.3 Targets 

The selected q uasars were known to be polarized in a broad 
optical band Ip tockma n et alJll984l ; iBerriman et al.lll990t 
IVisvanathan fc Willdll99Sft . We observed the ones with red- 
shifts larger than ~ 0.2 so that the Balmer edge in the ob- 
served frame is well away from the atmospheric cut off and 
well within the detector sensitivity range. Most of them are 
at redshift less than ~0.6 in order to cover the H/3 line. In 
the VLT observations, we also observed a few higher redshift 
objects to explore the Balmer continuum under the small 
blue bump. Four quasars in our sample (3C95, 4C09.72, B2 
1208+32, and 3C3 23.1) had already been observed spec- 
tropolarimetrically JSchmidt fc Smith! [2000 ) and had been 
found to have unpolarized broad lines with limited S/N. We 
have carried out long integrations on them to scrutinize the 
polarization of the broad emission lines and the Balmer edge 
waveleng ths. Ton 202 had been ob served spectropolarimet- 
rically bv lSchmidt fc Smith! i2000Tl and by us (Paper I), and 
was included in our sample for re-observation. In the VLT 
run, other quasars were observed with a short integration 
in the first night, and we tried to follow up in the second 
night the ones which seemed to show unpolarized broad lines 
and/or unpolarized small blue bump (see next section). 

2.4 Interstellar polarization 

Interstellar polarization (ISP) can affect our 
observations. For stars wit h E b-v ^ 0.1 mag, 
ISerkowski. Mathewson fc Ford! <ll975l) showed that 
Pmax<9EB-v where P max is the maximum polariza- 
tion in percent over the (optical) wavelengths. We list 



Table 4. Properties related to interstellar polarization. 



Name 


E 


B-V 


9E B _ V 




000034-1 58 





050 


0.45 


(107,-45) 


0090^4-094 





029 


0.26 


f1 58 55*1 


3C95 





.079 


0.71 


(205,-46) 


Q0350-073 





.079 


0.71 


(197,-43) 


Q0405-123 





.058 


0.52 


(205,-42) 


Q0414-060 





.043 


0.39 


(199,-37) 


Q1004+130 





.038 


0.34 


(225+49) 


B2 1208+32 





.017 


0.15 


(182+80) 


Ton 202 





.019 


0.17 


(36, +68) 


3C323.1 





.042 


0.38 


(34, +49) 


Q1912-550 





.059 


0.53 


(342,-25) 


Q21 15-305 





.120 


1.1 


( 16,-44) 


Q2251+113 





.086 


0.77 


( 83,-42) 


4C09.72 





.042 


0.38 


( 86,-46) 


Q2349-014 





.027 


0.24 


( 92,-60) 



E B -v is from NED (data from ISchlegel et al]ll99sh . 9E B -v is 
a maximum expected interstellar polarization along the line of 
sight to each quasar, but relatively uncertain for these low E B _y 
cases. (l,b) is Galactic coordinates. 

Eb-v and 9Eb-v as well as the Galactic coordinates (l,b) 
in Table [I] However, the limit of 9Eb-v are less certain 
for the lines of sight w ith Eb-v ^0.1 mag (see Figure 9 of 
ISerkowski et al.|ll975l i. A better constraint can be obtained 
from the polarization measurements of Galactic stars in 
each quasar field. We have check ed the polarization catalog 
of Galactic stars bv lHeilesl <l2000li . and have tried to observe 
a star for each quasar with a suspected ISP effect and/or 
for quasars observed with a long integration time. The stars 
were chosen to be nearby in the projected distance to each 
quasar but far enough in the real distance (essentially a 
faint early-type star or a faint giant) in order to sample 
the whole ISP in our Galaxy along the line of sight. We 
obtained the measurement for four quasars (4C09.72, 3C95, 
Q1912-550, 3C323.1), as summarized in Tabled 



3 RESULTS 

Our primary targets are the quasars with a small continuum- 
confined polarization, showing unpolarized broad emission 
lines. In these quasars, we expect in principle the small blue 
bump also to be unpolarized. Thus, these quasars would 
show a polarization which drops significantly at wavelengths 
shortward of 4000 A in the rest frame. This is a broad band 
property and relatively easy to recognize. We have measured 
this polarization drop as the ratio rp of the polarization at 
2891-3600A to that at 4000-4731A. Those are listed in Ta- 
ble El (pol arization in this table has been debiased following 
ISimmons fc Stewardll985h . 

We also have measured the polarization of the emission 
lines by fitting a power-law in the adjcent continuum region 
and subtracting off the continuum both in the total flux 
and polarized flux. More precisely, we first measured the 
PA of the continuum polarization, and rotated the Stokes 
parameters to that PA, and fit a power-law and subtracted it 
off. Then we measured the PA of the residual Stokes spectra, 
rotate the Stokes parameters to that PA, and measured the 
flux. Then we calculated the equivalent width (EW) of the 
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Table 5. The interstellar polarization measurements of Galactic stars around some of our quasars. Distances are 
from spectral type classifications, except for SA0245987 where it is from the measured trigonometric parallax. A max 
is the wavelength at which the polarization is maximum and P max is that maximum polarization at A ma x- 



Galactic star 


projected distance 


distance (kpc) 


-Pmax (%) 


Amax (A) 


PA(V) (°) 


HD23938 


1.2° NW from 3C95 


~ 0.57 


< 0.1 






AG+22 1545 


1.3° N from 3C323.1 


- 1.3 


0.718 ±0.005 


5570 ± 130 


84.3 ± 0.4 


SA0245987 


0.5° SW from Q1912-550 


0.24 ±0.08 


0.715 ±0.005 


5530 ± 110 


1.9 ±0.3 


Pn23sl-15 


0.4° NW from 4C09.72 


~ 3.3 


0.305 ± 0.007 


4970 ± 420 


104.7 ± 1.3 



Table 6. Broadband polarization properties. The measurements before any ISP correction are shown. For 
Q0003±158, 3C95, 3C323.1, Q1912-550, and 4C09.72, an additional row right below the row for each object presents 
the measurements after the ISP correction, with the first column showing the amount of ISP adopted. 



Name 


z 


2891-3600A 




4000-4731 A 


r P 






P (%) 


PA (°) 




P (%) 


PA (°) 




Q0003±158 


0.450 


0.72 ±0.02 


95.0 ±0.8 





.93 ± 0.02 


94.4 ±0.7 


0.77 ±0.03 


ISP 103°, 0.45% 




0.32 ±0.02 


83.8 ±1.7 





.52 ± 0.02 


87.0 ± 1.3 


0.61 ±0.05 


Q0205±024 


0.155 


0.38 ±0.09 


166.8 ±6.3 





.60 ±0.03 


3.0 ± 1.6 


0.62 ±0.15 


3C95 


0.616 


0.75 ±0.01 


176.0 ± 0.3 


1 


.17 ±0.02 


1.4 ±0.4 


0.64 ±0.01 


ISP 140°, 0.32% 




0.71 ±0.01 


8.6 ±0.3 


1 


.17 ±0.02 


8.9 ± 0.4 


0.61 ±0.01 


Q0350-073 


0.962 


1.01 ±0.04 


0.9 ±1.0 


1 


.03 ± 0.10 


170.0 ±2.7 


0.98 ±0.10 


Q0405-123 


0.573 


0.32 ±0.02 


129.9 ±1.8 





.42 ± 0.03 


138.9 ± 1.8 


0.77 ±0.07 


Q0414-060 


0.775 


0.80 ±0.05 


163.2 ± 1.8 





.90 ±0.10 


161.7 ±3.2 


0.88 ±0.12 


Q1004±130 


0.240 


1.43 ±0.02 


59.9 ±0.4 


1 


.17 ±0.01 


66.8 ± 0.3 


1.22 ±0.02 


B2 1208±32 


0.388 


1.01 ±0.01 


34.2 ±0.3 


1 


.41 ± 0.01 


28.9 ± 0.2 


0.72 ±0.01 


Ton202 


0.366 


1.25 ±0.01 


66.9 ±0.3 


1 


.87 ±0.02 


68.9 ± 0.2 


0.67 ±0.01 


Ton202 (May2002) 


0.366 


1.19 ±0.02 


179.5 ±0.4* 


2 


.11 ± 0.01 


0.2 ± 0.2* 


0.56 ±0.01 


3C323.1 


0.264 


0.78 ±0.01 


21.1 ±0.4 


1 


.37 ±0.01 


15.8 ± 0.2 


0.57 ±0.01 


ISP 90°, 0.90% 




1.51 ±0.01 


10.2 ±0.2 


2 


.19 ± 0.01 


9.6 ± 0.1 


0.69 ±0.01 


Q1912-550 


0.402 


1.14 ±0.02 


6.3 ±0.5 


1 


.75 ± 0.02 


4.5 ± 0.4 


0.65 ±0.01 


ISP 10°, 0.70% 




0.49 ±0.02 


1.1 ± 1.3 


1 


.09 ± 0.02 


1.0 ±0.6 


0.45 ±0.02 


Q2115-305 


0.979 


1.32 ±0.03 


62.0 ±0.7 





.78 ± 0.09 


55.3 ± 3.3 


1.69 ±0.20 


Q2251+113 


0.326 


0.83 ±0.07 


49.7 ±2.5 





.98 ± 0.06 


51.0 ± 1.7 


0.85 ±0.09 


4C09.72 


0.433 


1.02 ±0.01 


113.0 ±0.3 


1 


.33 ± 0.01 


119.2 ± 0.3 


0.77 ±0.01 


ISP 100°, 0.70% 




0.50 ±0.01 


131.4 ±0.6 





.93 ± 0.01 


131.4 ±0.4 


0.54 ±0.01 


Q2349-014 


0.174 


0.52 ±0.10 


144.5 ±5.2 





.89 ± 0.05 


154.7 ± 1.5 


0.58 ±0.11 



*The absolute PA for the May 2002 data of Ton 202 was not determined and only the relative PA was measured. 



polarized line flux using the fitted continuum (although the 
PA could be different) . The results are shown in Table |7| 
and Table |S] (quoted polarization in these two tables is the 
normalized, rotated Stokes q with the reference axis along 
the quoted PA, and thus has not been debiased). 

For the H/3 line, we measured the polarization and EW 
of the H/3 line at 4731 — 4934A, i.e. the wavelength region of 
[OIII]A4959,5007 lines were excluded to avoid the effect of 
these lines on the H/3 measurements. The EW thus derived 
is smaller than the real EW. This is not a problem for our 
primary purpose here, which is to determine whether or not 
the line is detected in the polarized flux and to determine 
its polarization if the line is detected. The EW of the line in 
total flux was measured using the same wavelength window, 
to compare directly with the EW in the polarized flux. The 
measurements for [OIII]A5007 line and MgIIA2798 line were 
done at 4982-5032A and 2705-289lA, respectively. 

Note that the absorption features at the observed wave- 
lenths of ~ 6900A and ~ 7600A seen in the total flux and 
polarized flux spectra are telluric. These regions were ex- 
cluded from the measurements (thus, also due to this exclu- 



sion, the EWs measured here can be smaller than the true 
EWs). 

In the following, we first describe the results for the 
quasars where the ratio rp is measured to be significantly 
less than unity with ~ 3a. The majority of the objects in 
our sample show this significant drop in P: 4C09.72, 3C95, 
3C323.1, B2 1208±32, Ton202, Q1912-550, Q0003+158, 
Q0405-123, Q2349-014. Then we briefly describe the results 
for the other six quasars: for these, we did not detect a sig- 
nificant P decrease shortward of 4000 A, although the con- 
straint is not tight in some of the objects. The polarizations 
of Galactic stars around each target other than those lis ted 
in Table are based on the compilation by Hciles (2000). 

All flux has been corrected for the Galactic reddening 
us ing the values of Eb-v in Ta bleland the reddening curve 
of lCardelli. Clayton, fc Mathisl Jl989D with R v = 3.1. 

3.1 4C09.72 

FigQ shows the observed polarization for 4C09.72. The up- 
per part with four panels shows from top to bottom (1) 
Stokes / (-Fa) j (2) normalized q' and u (thick and thin 
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Table 7. H/3 and [OIII] line properties. See the caption of Table for the ISP correction adopted. 



Name continuum region H/3 line region (4731-4934A) OIII5007 line region (4982-5032A) 

(4537-4731, 5032-5185A) polarization EW (A) polarization EW (A) 

pl%) pah p~(%) pA~n q 7 i pt%) pxn o 7 r 



Q0003+158 


0.92 ± 0.04 


95.2 ± 1.3 


0.53 ±0. 


.13 


94. ± 7. 


54. ± 13. 


92. 


0.40 ±0.23 


79. ± 17. 


14. ± 8. 


31. 


ISP 103°, 0.45% 


0.52 ±0.04 


88.7 ±2.2 


0.18 ±0. 


.13 


69. ± 20. 


33. ± 22. 




0.34 ±0.23 


42. ± 19. 


21. ± 14. 




Q0205+024 


0.60 ± 0.05 


4.1 ± 2.5 


0.75 ± 0. 


.37 


132. ± 14. 


51. ± 25. 


12. 


0.56 ±0.35 


46. ± 18. 


23. ± 14. 


25. 


3C95 


1.16 ±0.03 


2.0 ±0.8 


0.11 ±0. 


.17 


170. ± 42. 


5. ±7. 


50. 


0.48 ±0.24 


54. ± 15. 


8. ±4. 


20. 


ISP 140°, 0.32% 


1.17 ±0.03 


9.1 ± 0.8 


0.25 ±0. 


.16 


40. ± 19. 


11. ± 7. 




0.75 ±0.24 


54. ± 10. 


13. ±4. 




Q0405-123 


0.45 ±0.05 


138.3 ± 3.3 


1.12 ±0. 


.63 


96. ± 16. 


45. ±25. 


18. 


0.59 ±0.48 


55. ± 23. 


22. ± 18. 


18. 


Q1004+130 


1.11 ±0.02 


68.1 ± 0.4 


0.58 ±0. 


.26 


134. ± 13. 


9. ±4. 


17. 


0.56 ±0.32 


71. ± 16. 


4. ±2. 


8. 


B2 1208+32 


1.29 ±0.02 


28.1 ± 0.5 


0.29 ±0. 


.08 


171. ±8. 


17. ± 5. 


72. 


0.29 ±0.12 


36. ± 12. 


7. ±3. 


28. 


Ton202 


1.79 ±0.03 


69.5 ± 0.5 


0.33 ±0 


.11 


56. ± 12. 


11. ± 5. 


55. 


0.50 ±0.12 


79. ± 7. 


13. ± 3. 


40. 


Ton202 (May2002) 


2.07 ±0.02 


0.3 ±0.3* 


0.18 ±0 


.09 


3. ± 14.* 


6. ± 3. 


73. 


0.46 ±0.08 


8. ± 5.* 


11. ± 2. 


45. 


3C323.1 


1.42 ±0.02 


15.4 ±0.3 


0.91 ±0 


.08 


95. ± 3. 


33. ±3. 


51. 


1.19 ±0.13 


86. ±3. 


16. ± 2. 


18. 


ISP 90°, 0.90% 


2.24 ±0.02 


9.5 ±0.2 


0.19 ±0 


.08 


134. ± 13. 


4. ± 2. 




0.33 ±0.13 


78. ± 11. 


3. ± 1. 




Q1912-550 


1.89 ±0.04 


5.1 ± 0.6 


1.18 ±0 


.14 


175. ± 3. 


46. ± 5. 


73. 


1.15 ±0.16 


176. ±4. 


24. ± 3. 


39. 


ISP 10°, 0.70% 


1.23 ±0.04 


2.4 ± 0.9 


0.70 ±0 


.14 


160. ±6. 


41. ±8. 




0.65 ±0.16 


161. ± 7. 


21. ± 5. 




Q2251+113 


1.16 ±0.09 


50.1 ± 2.3 


1.75 ±0 


.39 


116. ±6. 


108. ± 24. 


72. 


1.13 ±0.63 


141. ± 17. 


21. ± 12. 


23. 


4C09.72 


1.32 ±0.02 


120.8 ±0.4 


0.56 ±0 


.16 


101. ±8. 


14. ± 4. 


33. 


0.65 ±0.24 


101. ± 10. 


7. ± 3. 


13. 


ISP 100°, 0.70% 


0.96 ±0.02 


132.5 ±0.6 


0.04 ±0 


.16 


1. ±90. 


1. ±6. 




0.02 ±0.24 


98. ± 90. 


0. ±3. 




Q2349-014 


0.91 ±0.07 


159.7 ±2.3 


0.80 ±0 


.43 


102. ± 15. 


45. ± 24. 


51. 


0.28 ±0.58 


42. ± 59. 


7. ± 14. 


21. 



*The absolute PA for the May 2002 data of Ton 202 was not determined and only the relative PA was measured. 



Table 8. Mgll line properties. Continuum region is 2611-2705A and 2891-2985A. Line region is 2705-289lA. See 
the caption of Table l6l for the ISP correction adopted. 



Name 


continuum region 


line region 


EW (A) 




P (%) 


PA (°) 


P (%) 


PA (°) 


Q' 


/ 


Q0003±158 


0.60 ± 0.06 


98.5 ±2.7 


0.65 ±0.26 


104. ± 11. 


41. ± 16. 


37. 


ISP 103°, 0.45% 


0.20 ±0.06 


89.2 ±8.1 


0.27 ±0.26 


100. ± 27. 


57. ± 54. 




3C95 


0.69 ± 0.02 


170.2 ± 0.8 


0.05 ±0.13 


82. ± 78. 


2. ±6. 


33. 


ISP 140°, 0.32% 


0.60 ±0.02 


3.4 ±0.9 


0.35 ±0.12 


55. ± 10. 


19. ±7. 




Q0350-073 


0.97 ±0.06 


178.4 ± 1.8 


1.03 ±0.38 


33. ± 11. 


33. ± 12. 


32. 


Q0414-060 


0.70 ±0.10 


159.8 ± 3.9 


0.58 ±0.43 


134. ± 22. 


36. ± 26. 


43. 


B2 1208±32 


1.00 ±0.04 


34.8 ± 1.1 


0.61 ±0.16 


136. ±8. 


23. ±6. 


38. 


Ton202 


1.10 ± 0.04 


66.4 ±1.1 


0.20 ±0.13 


165. ± 19. 


11. ±7. 


60. 


Ton202 (May2002) 


0.80 ±0.09 


175.0 ±3.3* 


0.45 ±0.16 


11. ± 10.* 


44. ± 16. 


78. 


3C323.1 


0.82 ± 0.05 


18.3 ± 1.8 


1.05 ±0.26 


79. ± 7. 


46. ± 11. 


37. 


ISP 90°, 0.90% 


1.52 ± 0.05 


9.4 ± 1.0 


0.50 ±0.26 


61. ± 15. 


12. ±6. 




Q1912-550 


0.98 ± 0.08 


6.2 ±2.4 


1.20 ±0.22 


5. ±5. 


73. ± 13. 


55. 


ISP 10°, 0.70% 


0.36 ± 0.08 


179.3 ± 6.5 


0.56 ±0.22 


177. ± 12. 


97. ±39. 




Q2115-305 


1.52 ± 0.05 


62.5 ± 1.0 


1.92 ±0.39 


72. ±6. 


34. ± 7. 


27. 


Q2251±113 


1.08 ± 0.47 


69.5 ± 12.3 


6.30 ± 1.71 


19. ±8. 


188. ±51. 


35. 


4C09.72 


0.94 ±0.03 


111.3 ± 1.0 


0.51 ±0.20 


105. ± 11. 


13. ±5. 


25. 


ISP 100°, 0.70% 


0.40 ±0.03 


132.4 ±2.3 


0.21 ±0.20 


175. ± 27. 


13. ± 12. 





*The absolute PA for the May 2002 data of Ton 202 was not determined and only the relative PA was measured. 



line, respectively), (3) unnormalized Q' and U' (q' x F\ and 
u x F\) with / (dotted line) scaled to match Q' at the red 
side, (4) position angle (PA) spectrum. The Stokes param- 
eters are rotated to have their reference axis along the PA 
of the polarization calculated at 4000-4731 A (Table HO this 
PA is indicated by the dotted line in the panel for the PA 
spectrum in FigHl. thus designated as q , etc. rather than 
q, etc. Therefore, q' is essentially the polarization P, and 
q' x F\ is essentially the polarized flux P x F\. We summa- 
rize the symbols used in this paper in Tabled The separate 
bottom part shows the polarized flux spectrum in vF v form, 
i.e. q' x vF v , against A with both axes in log scale (thus note 



that x-axes are different in these two parts). The total flux 
in vF v is also shown (dotted line). The rest wavelengths are 
shown at the top of each part. 



3.1.1 ISP correction 1 



1 The reader who wishes to skip to our best-corrected data can 
proceed to section [3.1.21 and Fig|2l ( and skip FiglSlandlH. Like- 
wise, in the following sections for each object, the ISP-related 
discussion is essentially separated in the first subsection when 
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Figure 1. VLT/FORS1 spcctropolarimetric measurements for 
4C09.72. From top to bottom, Stokes I, normalized q' and u' , 
unnormalizcd Q' and U' (the referece axis is along the polariza- 
tion PA at 4000-4731 A) with scaled I (dotted line), and PA. The 
fluxes are in units of 10 -18 erg cm -2 sec~ x A -1 . The separate 
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Figure 2. The ISP-corrected polarization of 4C09.72, assuming 
ISP PA of 100°. 
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0.000 
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-0.010 



-0.015 
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-0.015 



0.005 



Figure 3. The data on the plane of normalized Stokes parameters 
q and u with error bars. The numbers beside each data point are 
wavelengths in units of lOOA. The polarizations of H/3 and Mgll 
lines are also shown and labeled as /3 and Mg, respectively. The 
direction of the sky PA 0°, 45°, 90°, and 135° are also shown. 



One distinct observed feature in this object is that the PA 
rotates shortward of 4000A in the rest frame, probably flat- 
tens again shortward of the rest 3500A. We argue that this 
is most likely due to a contamination by ISP in our Galaxy 
and that the data can be corrected for this contamination 
to the spectra shown in Fig|2] 

In Fig|3] we put the same data shown in FigQon the 
q-u plane (normalized Stokes parameters; binned by 100 pix- 
els or 260 A). The size of the crosses shows the error. The 
labels are the central rest wavelengths in each bin in units 
of 100A. We see that the data points approximately align 
linearly, but the line of the alignment does not pass through 
the orgin (q,u) = (0,0). This is what we expect if the ob- 
served polarization consists of an ISP component and an 
intrinsic polarization of the object at a fixed PA. In general, 
if the polarization of an object changes in magnitude over 
the wavelengths but stays constant in PA, then the data 
points on the q-u plane align linearly with the alignment 
direction pointing toward the origin. On the other hand, a 
small contamination from an ISP, to a first order, will shift 
all the data points by a constant vector, since the ISP is in 
principle fixed at a certain position angle (PA) on the sky 
with its magnitud e not changing s ignificantly over the op- 
tical wavelengths dSerkowski et alJll975l) . The data points 
will stay linearly aligned, but the extrapolation of the align- 
ment direction will not pass through the origin in this case. 
This seems to be what we see in Fig[3] 

If the PA of the polarization of the emission lines is 
the same as the continuum or if the emission lines are in- 
trinsically unpolarized, the data point for the emission lines 
should also align along the other data. This is also the case, 



substantially addressed, so that the reader who wishes to skip it 
can proceed directly to the second subsection. 




4000 5000 6000 _ 7000 
Observed Wavelength (A) 

Figure 4. The data for 4C09.72 before any ISP correction. The 
upper panel shows the normalized Stokes parameter q with the 
reference axis at PA = 132° — 45°, i.e. perpendicular (on the 
q-u plane) to the PA of the quasar's intrinsic polarization. The 
smooth curve is the best fit Serkowski curve. The lower panel 
shows the unnormalizcd Stokes parameter Q with the reference 
axis at PA = 100° + 45°, i.e. perpendicular to the PA of the 
inferred ISP. 



as we show in Fig[3] the polarization of the H/3 and Mgll 
lines from Table |7| and |H| 

If this two-component interpretation is correct, the 
alignment direction of the data points on the q-u plane gives 
us the intrinsic PA of the polarization of the object. The ISP 
can roughly be inferred from the displacement of the data 
points: at least the displacement gives a constraint on the 
possible PA range of the ISP. However, the determination of 
an exact ISP vector requires an assumption of a PA of the 
ISP (or a magnitude of the ISP, where two possible solutions 
for the PA would exist). The polarization of Galactic stars 
in the same field can be used to check if it is consistent with 
the displacement of the data points seen on the q-u plane, 
and then to assume the PA of the ISP. 

The ISP probed by Galactic stars around 4C09.72 is 
consistent with this interpretation. Firstly, our observation 
of a star which is close to 4C09.72 in projection (see Ta- 
ble shows an ISP of 0.3% at PA 105°. This is overall a 
reasonable representation for the displacement of the origin 
of the quasar's intrinsic polarization at least in terms of its 
direction (see Fig|3J. Secondly, the observations of several 
stars around 4C09.72 in the literature compiled by iHeilesI 
( 2000) show a rather clear tendency for the PA of the ISP 
to be at 100°- 110° with P = 0.2 - 0.3% (for the stars of 
estimated distance of 100-300pc). Therefore, it is likely that 
the observation is simply affected by an ISP and the intrin- 
sic polarization of 4C09.72 is (at least approximately) at a 
fixed PA. 

For a given assumed PA of the ISP with a Serkowski 
curve (we also assume the empi rical relation of the constant 
in the Serkowski-law with Am a ^: lwiiittet et alll992T) . we can 
determine the best fit magnitude of the ISP which produces 
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the most flat PA spectrum for the intrinsic polarization of 
the object. Fig[5] shows the result of this correction. If the 
PA of the ISP is assumed to be at 100° (see below for the 
case of a larger PA), the best fit for the 4C09.72 data results 
in having a maximum polarization of P max = 0.70 ± 0.03% 
at A max = 4130 ± 540A for the ISP and 6» PA = 131.9 ± 1.7° 
for the intrinsic polarization of 4C09.72 (this PA was al- 
ready apparent from FigEJ ■ The Stokes parameters in Fig|5| 
have been rotated to have their reference axis to be at this 
intrinsic PA for the quasar. 

We can check if the ISP component in the observed po- 
larization really resembles a Serkowski curve, by looking at 
the spectrum of the Stokes parameter perpendicular in the 
q-u plane (or 45° in the sky plane) to the intrinsic polar- 
ization of the quasar, since this Stokes parameter spectrum 
should not contain the polarization of the quasar. This is 
shown in the upper panel of Fig[I] the spectrum is nicely 
fit by a Serkowski curve (this is of course essentially equiva- 
lent to checking whether the ISP-corrected stokes u, shown 
in Fig[5] is flat). Similarly, we can check the shape of the 
ISP-corrected polarized flux of the object by looking at the 
Stokes parameter perpendicular in the q-u plane to the PA of 
the ISP as shown in the lower panel of Fig|l] this essentially 
shows the same spectrum as the qF\ spectra in Fig|5| 

We note that the ISP magnitude derived above looks 
large compared with those for the nearby star we observed 
or several other stars in the literature. The estimated Eb-v 
for the line of sight to 4C09.72 is 0.042 mag, and thus we 
might expect that the ISP for 4C09.72 would be less than 
9E B -v ~ 0.4% (see Table 0}. However, the flatness of the 
PA spectrum after the ISP subtraction (Fig[5J makes this 
correction quite believable. We note again that the limit of 
9Eb-v are less cert ain for the lines of sig ht with smaller 
Eb-v in Figure 9 of ISerkowski et alJ il975l> . 

The EW of the emission lines in the corrected polarized 
flux is consistent with zero and its upper limit for the H/3 
line is much less than the EW of the line in the total flux 
(Tables [7| and |HJ. This is what we expect since the assumed 
PA of the ISP is essentially the same as the observed PA 
of the line polarization, and the data points for the lines in 
the q-u plane are essentially aligned with other data points. 
Conversely, our interpretation is well supported by this line 
polarization. 

The uncertainty in the PA of the ISP leads to an un- 
certainty in the corrected polarized flux spectrum. However, 
the PA of the ISP is not likely to be much smaller than 100° 
based on the ISP map around the quasar. Conversely, if we 
increase the PA assumed, the features seen in the corrected 
polarized flux which are described below become stronger 
in this particular case (this makes the magnitude of the in- 
ferred ISP even larger as we can see from FiggJ. Thus, the 
overall spectral features do not depend on the ISP correc- 
tion at least qualitatively. The uncertainty in the PA of the 
ISP may be minimized by the future observations of a few 
more close- by distant stars around 4C09.72. 



3.1.2 Corrected spectra 

FigE?l shows the spectra for 4C09.72 which have been cor- 
rected for the contamination from the ISP in our Galaxy. 
Note that the S/N of the polarization spectra has been ef- 



fectively degraded due to the reduction of the P magnitude 
from the subtraction of the ISP. 

The corrected polarized flux (unnormalized Stokes Q' , 
or q' x Fx) shows spectral features that look similar to those 
seen in Ton 202 (Paper I). Firstly, there are essentially no 
emission lines, which are quantified in Table [7] and |H| Sec- 
ondly, in the Balmer edge region, there is a slope down-turn 
(or local maximum) at about 4000A, and there also seems 
to be a slope up-turn (or local minimum) at ~ 3600A. In 
addition, there seems to be a local peak feature at ~3050A. 
These are all discussed in section 0] 



3.2 3C95 

3.2.1 ISP correction 1 

The observed polarization for 3C95 is shown in Fig|5] We 
have combined the observations with an order sorting fil- 
ter GG435 (thus A bs > 4600A) and the ones without any 
filter (thus A obs < 6400A). There is again a PA rotation 
shortward of rest 4000 A. This is probably due to a small 
ISP contamination and we show the data corrected for the 
ISP contamination in Fig|S] This inference is again based 
on the q-u plane plot where the data points align at least 
approximately linearly as shown in Fig|7| (however see be- 
low) . The polarization measurements of a few Galactic stars 
having estimated distance of 100-600pc sug- 
gests an overall tendency of the ISP at PA 130-150°with 
P — 0.1 — 0.4%, though our observation of a neighboring 
star (Table|KJ shows no polarization larger than 0.1%. This 
direction of the PA seems to be a reasonable representation 
for the shift of the intrinsic polarization of 3C95 on the q-u 
plane (Fig0. 

Thus the argument for ISP contamination looks rea- 
sonable. However, there also seems to be some systematic 
deviation from the alignment on the q-u plane at the blue 
end of the spectrum, shortward of the rest ~3000A. There- 
fore, we tried to correct for the ISP using the data long- 
ward of this wavelength. The same correction procedure 
as for 4C09.72, assuming the PA of the ISP as 140°, pro- 
duced the results shown in Fig|S] The fit resulted in having 
Pmax = 0.32 ± 0.03% and the PA of the quasar's intrin- 
sic polarization at 9.7 ± 0.9° (A max was simply assumed to 
be 5500A). The correction looks overall reasonable except 
for the blue end (the ISP correction using the whole wave- 
length results in having some wave-like rotation in the PA 
spectrum) . 

There is an uncertainty in the polarized flux shape due 
to the uncertain ISP PA, but this is rather small since the 
magnitude of the ISP is small compared to the observed 
polarization. If the PA of the ISP is assumed to be smaller 
than 140° , this would give a larger ISP magnitude as can be 
seen in Fig|7| which is unlikely as we argued above. For an 
ISP correction with a larger PA, the spectral features in the 
polarized flux even become slightly stronger. 

The H/3 and Mgll lines in the spectrum before any ISP 
correction are unpolarized with a ~ 0.15% (see Tabled and 
[SJ. This means that the ISP contamination is small (0.3- 
0.4% at most) and the lines are intrinsically unpolarized, 
unless the ISP is coincidentally canceling out any intrin- 
sic polarization component. (The above ISP correction 'pro- 
duces' a little line polarization, but not significantly). 
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Figure 5. The same as Fig0 but for 3C95 



Figure 6. The ISP-corrected polarization of 3C95, assuming ISP 
PA of 140° . 
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Figure 7. The same as Fig|3] but for 3C95. 



3.2.2 Corrected spectra 

Figlol shows the data corrected for the ISP contamination. 
The polarized flux spectrum again shows a spectral feature 
at the Balmer edge region similar to that of Ton 202 in Paper 
I and 4C09.72 in the previous section: a slope down-turn at 
around rest 4000A and an up-turn at around 3600A. As we 
argued above, the emission lines probably are intrinsically 
unpolarized, even though the [OIII] line and Mgll lines look 
marginally polarized in the corrected polarized flux spec- 
trum. 



3.3 3C323.1 

3.3.1 ISP correction 1 

The data for this object, shown in Fig|5] are also likely to 
have been affected by ISP: the data corrected for the ISP 
are shown in Fig[^] 

The original measurement in Fig|H] shows a PA rota- 
tion shortward of rest 4000A and also at emission line wave- 
lengths, but on the q-u plane, a linear alignment of the data 
points is seen (Fig llUH . The direction of the shift of the data 
points is consistent with the observations of Galactic stars 
in the same field. Our observation of a nearby star shows 
P = 0.7% at PA 84° (see Table d, and other observations 
of three Galactic stars within 3° (Hcilcs 200j3) having esti- 
mated distance larger than 200pc also show a similar PA 
(85-95°), though with a smaller P (0.2-0.3%). Therefore we 
have corrected the data assuming the PA of the ISP at 90° , 
and Figgis the result of this correction. Before this ISP cor- 
rection, there were some apparent features at the emission 
line wavelengths in the polarized flux (FiglHj, but after the 
correction, these are essentially gone. This is expected since 
the polarization measured for these lines before the correc- 
tion is roughly consistent with the subtracted ISP fFig llOl 
Tables 00. 

The uncertainty in the PA of the ISP leads to the un- 
certainty in the polarized flux. An ISP correction with a 
smaller PA results in a stronger spectral feature in the po- 
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Figure 10. The same as Fig|3] but for 3C323.1. 



larized flux. With a larger PA, the feature becomes weaker, 
but this also results in having emission lines in the polarized 
flux. The ISP magnitude also increases rapidly with a larger 
assumed PA. 



3.3.2 Corrected spectra 

The data which are best-corrected for the ISP are shown in 
FiglHl The emission lines look unpolarized (see also Table |7| 
and|8|) after the ISP correction. At the Balmer edge region, 
there remains the now-familiar intrinsic features left: the 
slope down-turn at ~ 4000A, and an up-turn at ~ 3600A, 
though the feature looks weaker than for the previous ob- 
jects. 



3.4 B2 1208+32 

3.4.1 ISP effect 1 

In this object, a small PA rotation shortward of 4000A is 
again seen fFig ll IP and the data points on the q-u plane 
also roughly align linearly (Fig ll2L though not too cleanly. 
We do not have our own observation of a Galactic star in the 
field, but the data in the literature indicate that stars with 
estimated distances of 400-500pc show a small polarization 
(<0.25%) at PA 0-25°. 

If there is any effect of ISP around at this PA, it is ex- 
pected to be small, judging from the location of the data 
points on the q-u plane. However, we note that the contin- 
uum PA rotation direction and the line polarization direc- 
tion are not in the same sense, which suggests that one of 
them is intrinsic. In either case, the ISP magnitude is prob- 
ably small, so its effect should also be small. 

3.4-2 Uncorrected spectra 

Fie llll shows the measurements without any ISP correction, 
and the intrinsic polarization is expected to be qualitatively 
the same, since the ISP effect is thought to be small. The 
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Figure 12. The same as Fig|3] but for B2 1208+32. 



measured small line polarization (polarized at different PA 
from the continuum) could either be intrinsic or entirely 
due to an ISP. The polarized flux shows a similar, now- 
familiar feature in the Balmer edge region, and the up-turn 
at ~ 3600A in this object looks more pronounced. 



3.5 Ton202 

This object was found to show a Balmer edge feature in the 
polarized flux, when observed in May 2002 (Paper I). How- 
ever, we do not find such a dramatic feature in our data 
taken in May 2003, as we show in Fig |13l The polarized 
flux looks quite featureless at ~ 4000A. (At the red end 
wavelengths longer than ~6400A, the polarized flux shape 
become convex, which could be due to second-order light.) 
This difference in two epochs seems to be due to a real vari- 
ation in polarization. 

The upper panel of Fig |14l shows the difference in po- 
larization in these two epochs (more accurately, difference 
in normalized Stokes q rotated to the object's polarization 
PA at each epoch; note that the sky PA of the May 2002 
observation could not be determined due to some calibra- 
tion problem; see Paper I). Thus the polarization seems to 
have decreased longward of ~ 3700A, while the continuum 
flux shape stayed roughly the same as shown in the lower 
panel of Fig |14l (absolute flux levels have uncertainties from 
slit loss and seeing size difference, which are rather hard to 
quantify) . 

Polarization variability on a timescale of one year is not 
unexpected, since the polarization is thought to originate 
from a very small spatial scale, smaller than the size scale 
of the broad line region, as we argue in Paper I and also in 
the next section. In fact, polarization variability of Ton 202 
an d some other quasars in o ur sample has been pointed out 
bv lSchmidt fc Smith! feOOOl) . Some broad-line ra dio galaxies 
are a lso known to show polarization variability jAntonuccil 
Il984h . Ton 202, as well as other objects, should be followed 
up to confirm and explore the variability. 

We note that the [OIII]A5007 line might be weakly po- 
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Figure 14. Comparison of the observations of Ton202 at two 
epochs. The upper panel show the polarization observed in May 
2003 minus that in May 2002. The lower panel show the total flux 
in May 2003 divided by that in May 2002. 

larized. The measurements in two epochs gave consistent 
results for the magnitude of polarization (Table 0. Assum- 
ing that this [OIII] line polarization did not vary, we would 
be able to constrain that the PA of the continuum polar- 
ization has changed only less than ~ 9°(1ct) based on the 
la uncertainty for the [OIII] line polarization PA (Table 0. 
There might be a very small rotation in PA across 4000A 
(FigESJ, which was not seen in the data of May 2002, but 
no useful ISP measurement is available at this time. 

3.6 Q1912-550 

3.6.1 ISP correction 1 

As shown in Fig ll5l the polarization clearly decreases short- 
ward of rest 4000 A, so the small blue bump is at least par- 
tially unpolarized. However, the H/3 and [OIII] lines are po- 
larized as seen in the q x Fx spectrum. The data are proba- 
bly affected by the ISP, and we show the ISP-corrected data 
in Fie llfcil The data points in the q-u plane shown in Fig |17l 
align roughly linearly, suggesting the PA of the intrinsic po- 
larization is roughly 0-10°. Our observation of a nearby star 
shows a rather strong ISP of ~0.7% at PA ~2° (Table 0. 
Therefore, the observation of the quasar is probably affected 
by the ISP which happens to be at a PA similar to that of 
the object polarization, causing very little PA rotation in 
the observed data. 

The Galactic latitude is not high for this object, though 
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Figure 15. The same as FigQ] but for Q1912-550. 



Figure 16. The ISP-corrected polarization of Q1912-550. 
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Figure 17. The same as Fig|3] but for Q1912-550. 



the estimated Eb-v is not much larger than other objects 
(see Table0J ■ The existing polarizat ion measurem ents of the 
nearby Galactic stars (within ~ 6°; Heilcs 2000) show that, 
while two stars have low P ( < 0.15%) at PA 70-80° (with the 
estimated distance of 90-250pc), three stars have a rather 
high P of 0.4-0.6% at PA 5-15° (distance 170-440pc). This is 
consistent with our observation of one Galactic star. 

Fig ll6l is the result after the subtraction of ISP assum- 
ing PA=10°(with the same A max as the observed star in 
Table In this case, the best fit ISP magnitude was of 
fmax = 0.7%. With this correction, the narrow lines are 
rather well removed from the polarized flux spectrum, but 
there appear to be broad components (these were seen be- 
fore the correction) , which is quite different from the cases of 
4C09.72 and 3C95 or 3C323.1. The ISP magnitude could be 
slightly larger, if its PA is smaller than 10° (though the mag- 
nitude would not be constrained well when it is too close to 
the PA of the intrinsic quasar polarization) ; the broad lines 
still remain even with this larger ISP subtraction. In fact, 
this seems to be due to a qualitative difference in the original 
polarization spectrum, rather than due to undersubtraction 
of the ISP — unlike those seen in 4C09.72 (and in 3C95, 
3C323.1, and B2 1208+32), the polarization does not seem 
to decrease at the wavelengths of broad line wings. There- 
fore, the broad lines seem to be slightly polarized. 

Since the ISP PA seems to be very close to the PA of 
the quasar polarization, this leaves a corresponding rather 
large uncertainty in the broad-band shape of the polarized 
flux. 



3.6.2 Corrected spectra 

Fig ED shows the ISP-corrected spectra. In this object, the 
broad line components seem to be slightly polarized. One 
would then expect the small blue bump also to be somewhat 
polarized. Therefore, the decline of the intrinsic polarized 
flux continuum shortward of 4000A could be sharper than 
seen in the corrected polarized flux. However, we note that 
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Figure 18. The same as Fig[Q but for Q0003+158. 
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Figure 19. The ISP-corrected polarization of Q0003+158. 
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Figure 20. The same as Fig|3] but for Q0003+158. 



the uncertainty in the broad-band shape of the corrected 
polarized flux for this object is rather large. 



3.7 Q0003+158 

3.7.1 ISP correction 1 

This object fFig |18H also has an at least partially unpo- 
larized small blue bump, since P decreases shortward of 
4000A, but the broad lines look partially polarized. There 
is no clear evidence of PA rotations. Based on the polar- 
ization measurements of two Galactic stars ar ound this ob - 
ject (0.06/0.41% for 140/280pc at 115/118° : Hfeilel l2000fl . 
there could be a small ISP at a PA very similar to that 
observed for Q0003+158. Therefore, at least partially, the 
polarized lines could be due to ISP. In Fig |19l we illustrate 
a possible "nominally" maximum effect of an ISP correc- 
tion: namely we implemented the same fitting procedure as- 
suming P m ax = 0.45% which is equal to 9Eb-v, instead of 
assuming the PA of the ISP. In these corrected data, the 
polarized flux spectrum shows almost no emission lines (i.e. 
the emission line polarization seen in the uncorrected data 
is consistent with the assumed ISP; see Table |7| and [S] as well 
as the data on the q-u plane shown in Fig |20t , 



3.7.2 Corrected Spectra 

Fig |19l shows the ISP-corrected spectra, but note that we 
do not have a firm basis for the ISP estimate for this ob- 
ject. In these corrected spectra, almost no emission lines are 
seen in the polarized flux, and at the Balmer edge region, 
there seems to be again a slope down-turn at rest ~4000A. 
However, we need more data to confirm these, and also need 
more data on Galactic stars around this object. 
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Figure 21. The same as FigQ] but for Q0405-123. 
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Figure 22. The same as Fig|3] but for Q0405-123. We calcu- 
lated and plotted the line polarization of H/3+[OIII] region (4731- 
5032A) instead of H/3 only, since excluding the telluric absorption 
wavelengths at around 7600A leaves a very small wavelength win- 
dow for H/3 line. 



3.8 Q0405-123 

3.8.1 ISP effect 1 

A PA rotation across 4000A is seen (Fig l2H which resembles 
that of 4C09.72 (and 3C95 or 3C323.1). This could be due 
to ISP since the data points in the q-u plane roughly align 
linearly, with the alignment line displaced from the origin 
(Fig |22L although no good polarization measurements for 
Galactic stars in the quasar field are available. The align- 
ment direction suggests the intrinsic polarization is at PA 
~ 160°. The line polarization is at least consistent with 
this interpretation, but does not provide significant limits 
(FigHU the polarization of the H/3 + [OIII]A5007 line mea- 
sured at 4731-5032A is plotted instead of H/3 only, since ex- 
cluding the telluric absorption wavelengths at around 7600A 
left a very small wavelength window for H/3). 



3.8.2 Uncorrected spectra 

The data without any ISP correction are shown in Fig |21l 
The overall polarization is very low (Fig |21l note the small 
P of ~ 0.5% even at the red side), but we marginally de- 
tected a polarization decrease shortward of 4000A (Tablc^. 
The above consideration of a possible ISP effect might sug- 
gest that the intrinsic polarization of this object is at PA 
~ 160°. The lines are consistent with being unpolarized, 
but the constraints are not tight. 

The object shows variability, as discussed by 
ISchmidt fc Smith! |2000l): when observed in 1990-1991 
bv~ fA^itonu^c^t^iIIl^96l) . the optical polarization was 
much lower than we observed (0.04 ± 0.02%). The mea- 
surement bv lSchmidt fc Smith! <2000Jl in 1995 and 1999 did 
not show significant optical polarization (P — 0.27 ± 0.19% 
and 0.17 ± 0.14%, respectively, though stil l consistent with 
our measurement. The measurements by IStockman et alj 
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(1984) in 1978 and 1980 are roughly consistent with ours, 
or P was slightly higher in one occasion in 1978. 

3.9 Q2349-014 

Although the formal polarization ratio rp is less than unity 
at 3 a (Table® , the polarization seems to behave differently 
from the other objects fFig |23H in that it looks constant over 
4000A, but rather abruptly decreases shortward of ~ 3500A. 
This interesting behavior should be confirmed and followed 
up with longer integrations. There could be a small effect 
from an ISP, since the polarization measurements of several 
Galactic stars show a rather clear PA alignment at 120-130° 
with P = 0.2 - 0.4%. 

3.10 Q0414-060 

We did not detect a significant P decrease shortward of 
4000 A (FigEU, and the limit on the Mgll line polarization 
is not significant either (Table®. No good polarization mea- 
surement for Galactic stars is found in the literature within 
5° from the quasar. 

3.11 Q1004+130 

3.11.1 ISP or second polarization component 1 

The polarization gradually increases toward shorter wave- 
lengths in this object, and the PA rotates also gradually 
(Tig |251 . The data points on the q-u plane align linearly 
(Fig |26L and its alignment probably suggests that it has 
one polarization component at PA ~ 45° and another com- 
ponent with P being roughly wavelength-independent and 
larger than ~ 0.8% based on the distance between the align- 
ment line and the origin (for this lower limit case, the PA 
of this component is at PA ~ 90°). The latter component 
might be an ISP: we do not have stars with large estimated 
distance within ~ 4° , but at radii 4 — 5° , several stars show 
a polarization roughly aligned at PA 75-85° , though the po- 
larization are all low ~ 0.1% (one has an estimated distance 
of lOOpc). We need a good observation of a Galactic star to 
check the ISP. 



3.11.2 Polarization and the BAL property 

In any case, the polarization property is quite different from 
the previous objects. The polarization seems to consist of 
two components: one of them is probably quite wavelength- 
independent, but the other has quite a blue polarized flux, 
which might be due to dust scattering. 

This quasar is one of the few Broad Absorption 
Line (BAL) quasars having both a well-defined radio 
jet axis and good spectropolarimetric data. The evidence 
for the BAL property of t his quasar was presented in 
I Wills. Brandt, fc Laorl jl999f) . The spectro polarimetry data 
of thi s quasar has been documented bv lAntonucci et alJ 
(1996), and our data are qualitatively consistent with their 
results in terms of the polarization magnitude and PA rota- 
tion, but their integrated polarization PA (though covering 
up to longer wavelengths) is a little larger than ours (by 
~ 10°), which might be due to variability. 

The overall alignment between the polarization PA 
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Figure 23. The same as Fig[Q but for Q2349-014. 
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Figure 24. The same as FigQ] but for Q0414-060. 



Figure 25. The same as Fig0 but for the Keck observation of 
Q1004+130. 
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Figure 26. The same as Fig0 but for Q1004+130. 



(67°) and radio jet axis (121°) is intermediate. However, 
the component at PA ~ 45° is roughly perpendicular to the 
radio axis, and the second component might be due to an 
ISP as described above. 

In the other two cases of BAL quasars with FRII 
jet axis and opti cal spectropolarim etric data (FIRST 
J101614.3+520916. iGregg et all l200Ct LBQS 1138-0126, 
iBrotherton et al.ll2002l) . the relations between the polariza- 
tion PA and jet axes are again intermediate and PA rotations 
are seen. For both cases, P is quite large and the ISP seems 
to be quite small based on very small Eb-v, suggesting that 
the PA rotation would originate from two intrinsic polariza- 
tion components rather than from an ISP contamination. 



3.12 Q2251+113 

The polarization looks constant over all the wavelengths, 
except at the emission lines (Fig |27t . The blue side of the 
H/9 line seems to be polarized at a different PA from the 
continuum. However, there might be some effect from an 
ISP, since the polarization measurements of several stars 
within 5° from the quasar show an ISP at a consistent PA 
of 100-105°with P = 0.2%. 



3.13 Q2115-305 

Polarization is rising towards the bluer side, and thus the po- 
larized flux is bluer than the total flux (FigUHJ. The Mgll 
line is polarized in the same way as the neighboring contin- 
uum (Table|HJ . This may well be a case of scattering by dust 
outside the BLR. Several stars around the quasar (within 4°) 
show a small polarization (P = 0.1 — 0.4%), and the PA ten- 
dency is not clear although two of them with large estimated 
distances (2.4kpc and 0.5kpc) have PAs rather close to the 
polarization PA of Q2115-305. 
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Figure 27. The same as Fig[Q but for Q2251+113. 
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Figure 28. The same as Fig[I] but for Q2115-305. The refer- 
ence axis of q' and «' is taken along the PA of the polarization 
integrated over the whole wavelength region, rather than at 4000- 
4731A. 



Figure 29. The same as FigEJ but for Q0350-073. The refer- 
ence axis of q' and u' is taken along the PA of the polarization 
integrated over the whole wavelength region, rather than at 4000- 
473lA. 
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Figure 30. The same as FigQ but for Q0205+024. 



3.14 Q0350-073 

The polarization looks almost constant, though slightly 
higher around rest 4000A, and the polarized flux shape un- 
der the small blue bump looks the same as the total flux 
(FiglH). There seems to be some PA rotation at the red 
end of the spectrum. Mgll line might be polarized at a dif- 
ferent PA from the continuum. Only small polarizations for 
Galactic stars (P < 0.2%) are found within 10° from the 
quasar, and no clear PA tendency is seen. 



3.15 Q0205+024 

This is a quasar which shows the spectrum of narrow-line 
Seyfert 1 galaxies. The polarization level is small (Tig l30H . 
It looks constant over 4000A and the ratio is consistent with 
unity but with a rather large uncertainty. Two stars (within 
5° from the quasar) show 0.2-0.3% polarization at 140-160°. 



4 DISCUSSION 

4.1 The Balmer edge feature 

4-1.1 Observed edge feature 

In Fig l31l we show the comparison of the q' 'Fx spectrum, i.e. 
essentially the polarized flux spectrum, for Ton 202 from 
Paper I and for 4C09.72, 3C95, 3C323.1 (after the ISP cor- 
rection) and B2 1208+32. The wavelengths are all in the rest 
frame of each quasar. The q'vF v spectra for all five objects 
are compared in Fig l32l Note that in the ISP corrections, 
the uncertainty in the PA of the ISP leads to some uncer- 
tainty in the ISP-corrected polarized flux shape, but it is 
qualitatively unchanged as we have described in the previ- 
ous section (though the shape in 3C323.1 is a little more 
ambiguous than in 4C09.72 and 3C95). 

The comparison shows that the slope of the polarized 
flux of all five quasars longward of ~4000A is essentially the 
same as that of the total flux continuum, but the slope has 
a down-turn at ~4000A (i.e. a local maximum at ~ 4000A: 
we note in one case, 3C323.1, this seems to be at a sub- 
stantially shorter wavelength than in others), around where 
the small blue bump starts in the total flux. In addition, 
another change (an up-turn) of the slope is seen at around 
3600A (i.e. the spectra have local minima at ~ 3600A). This 
slope up-turn looks more evident in B2 1208+32, while both 
of the slope changes look very weak in 3C323.1. The slope 
up-turn was already suggested in Paper I for Ton 202, but 
the comparison of these five quasars rather strengthens the 
suggestion, since the position of this slope up-turn looks 
roughly similar in the five polarized flux spectra. 

Based on all these, we identify this spectral feature, 
which is composed of these two slope chang Balmer 
edge absorption feature: (1) the slope change at ~ 4000A 
would be the start (long wavelength limit) of a broadened 
Balmer edge absorption feature (2) the change of the slope 
at ~ 3600A would be the position of the maximum opacity, 
with the flux shape at the shorter wavelength side being 
formed by the combination of the intrinsic continuum color 
and the v~ 3 decrease of the continuum opacity, although 
the slope looks still very red in the first three objects. The 
feature may be widespread, since the edge features in the 
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Figure 31. Comparison of the polarized flux spectrum, or q' F\ Figure 32. The same as Fig l31l but in vF v with both axes in log 

spectra, for Ton 202, 4C09.72, 3C95, B2 1208+32 and 3C323.1 scale, 
(solid lines). The spectrum of Ton 202 is from the data in Paper 
I. Flux has been corrected for Galactic reddening. For 4C09.72, 
3C95 and 3C323.1, the spectrum is after the ISP correction, i.e., 
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five objects look rather similar, and a few other objects show 
behavior which may also be similar, but which is less clear. 

In addition, there seem to be more spectral features in 
the shorter UV wavelengths. We will briefly describe those 
in section l4~2l 



4-1-2 Simplest Interpretations 

The polarized flux of all these quasars essentially show no 
emission lines, and thus the polarization is confined to the 
continuum. These quasars have lobe-dominant radio jets, 
and the direction of the optical continuum polarization ( El- 
vector direction) is in all cases parallel to the radio struc- 
ture, as we discuss in the next section. The unpolarized 
broad emission lines in our quasars lead to the inference that 
the polarization originates interior to the BLR. 2 The exact 
cause of the polarization is not well understood in these 
quasars, but we have discussed some simple possibilities in 
Paper I, all based on some form of electron scattering. Dust 
scattering is unlikely, since the region is thought to be inte- 
rior to the BLR and thus within the dust sublimation radius. 
The electron scattering might originate in the atmosphere of 
putative plane-parallel scattering-dominated disk. However, 
this case suffers from the wrong prediction of the polariza- 
tion direction, i.e. perpendicular to the symmetry axis of 
the disk which is thought to be along the jet direction (see 
Paper I for more discussion). Alternatively, the polarization 
is caused by an optically thin oblate pure electron scattering 
region surrounding the BBB emitter (with the BBB emis- 
sion being unpolarized, e.g. completely depolarized by Fara- 
day de polarization with magnetic fields in the disc atmo- 
spher e; lAgol fc Blaesl II99& lAeol. Blaes fc Ionescu-Zanettil 
1998). This produces the right polarization direction if the 
minor axis of the oblate shape is along the jet. In this latter 
case, the polarized flux is simply proportional to the BBB 
emission. In the case of electron scattering in the disk at- 
mosphere, the polarized flux is being produced as a part of 
the BBB emission and is thus intrinsic to the BBB emission. 
Therefore, in both of these cases, the Balmer edge feature 
seen in our quasars should be intrinsic to the BBB emission. 

To summarize, our most simple interpretation is that 
the observed Balmer edge absorption feature is intrinsic to 
the BBB emission. If this is correct, then the observed fea- 
ture directly indicates that the BBB emission is indeed ther- 
mal, and its emitter is optically thick. However, there is 
much to learn from the details of the polarized flux features, 
which tend to repeat from object to object, if we can identify 
them in theoretical models. We also discuss other possibili- 
ties below in section 14.1.41 



2 These properties are quite in contrast with Seyfert 2 galaxies 
where the broad emission lines and continuum are both polarized 
in the same way, and the pola rization is very high and perpen- 
dicular to the radio structure iAntonuccilll993l) . Based on these 
properties, the polarization in these type 2 objects is thought 
to be caused by polar scattering of the radiation from a hidden 
nucleus and BLR along the jet direction by matter in a region 
outside the BLR,. 



20 



^ 15 



5 




ML .. i .... i .... i .... i .... i ... I 

3000 3500 4000 4500 5000 
Rest Wavelength (A) 

Figure 33. As an illustration, a spectrum (thick smooth line) 
for a face-on disk model around a Kerr (a/M = 0.998) black 
hole of mass 8 X 10 9 Mq, accreting at 2 Mq yr — 1 is plotted 
over the May 2002 polarized flux spectrum of Ton202 (corrected 
for the Galactic reddening). The model spectrum has arbitrary 
normalization. 

4-1-3 Broadening and a disk atmosphere model 

The onset of the edge feature looks broadened over the whole 
range from ~ 360uA to ~ 4000A. This could be at least 
partly due to the high-order Balmer absorption lines and 
possibly metal lines as well. However, the observed feature 
suggests a broadening due to a high velocity dispersion in 
the gravitational potential well. Very naively, in terms of 
a simple black-body disk, the expected velocity dispersion 
looks roughly consistent with the observed extent of broad- 
ening. The radius r(4000A) emitting the radiation around 
4000A is given as 

r(4000A)/r 9 

where M, rj and ^Bdd are the mass of the black hole, ra- 
diative efficiency (L = r/Mc 2 with accretion rate M) and 
Eddington luminosity, respectively, and r g = GM/c 2 . The 
radius 350Y 9 corresponds to v ~ 0.05c for a Keplerian rota- 
tion, thus roughly consistent with the apparent broadening. 

Detailed disk model spectra (Hubeny et al. 2000) are 
roughly in agreement with these simple estimates, although 
considerable flux in the Balmer edge region can arise from 
smaller radii in some disk models. As mentioned in section^ 
accretion disk modeling of the Balmer edge region is com- 
plicated by the fact that the relevant annuli contain internal 
hydrogen ionization zones. These result in very steep, con- 
vectively unstable temperature gradients, as well as density 
inversions. This only adds to the considerable uncertainty in 
the vertical structure of such models. 

As an illustration, Fig l33l shows a face-on disk model 
around a Kerr (a/M = 0.998) black hole of mass 8 x 10 9 M , 
accreting at 2 Mq yr _1 . The model spectrum has arbitrary 
normalization, and is plotted on top of the Ton 202 Keck 
data from May 2002. Disks viewed at higher inclinations will 
have the edge smeared out further and shifted toward the 
blue. The model does not include any lines, particularly high 
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Figure 34. In the upper panel, the five polarized flux spectra 
shown in Fig l31l normalized by the mean at 2800-5000A, are over- 
plotted with a larger wavelength bin of 20A. In the lower panel 
the average of these five spectra with equal weighting is shown. 
Both panels are in F\. 



order Balmer lines and metal lines that may help shift the 
edge to longer wavelengths. Note that the model spectrum 
is the total flux spectrum, and some polarizing mechanism 
would still need to be invoked. (The intrinsic polarized flux 
produced by the disk model at nonzero inclinations has a 
much stronger edge feature, and is polarized in the wrong 
direction.) 

If the Balmer edge arises from a disk, an atmosphere 
contributing to the edge at some radius might change on a 
time scale comparable to the thermal time at the effective 
photosphere. This will be less than the local thermal time 
of the entire disk thickness at that radius, which is 

r~0.2yr(^) " (^) (^) . 

where a is the Shakura & Sunyaev (1973) viscosity param- 
eter. This might be consistent with the time scale between 
observations. However, the disk would have to do this in a 
coordinated fashion at all radii producing flux in the Balmer 
edge region, and this probably requires reprocessing of radi- 
ation from the inner disk. (Such reprocessing is already re- 
quired in disk models in order to explain the near simultane- 
ous optical/UV continuum variability observed in quasars.) 
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Figure 35. The same as Fig l34l but in vF v with both axes in log 
scale. 



4-1-4 Alternative possibilities 

As discussed in Paper I, we cannot rule out some alternative 
explanations. One possible case would be that the feature is 
due to foreground absorption. However, a simple foreground 
absorption is unlikely, due to the lack of a huge Lyman edge 
absorption in quasars. This is expected in a simple fore- 
ground absorption case, since the Lyman edge is a resonant 
feature, unlike the Balmer edge. In particular, for the five 
objects shown in Fig lbTI the ones with available far-UV spec- 
tra ( 4C09.72 and Ton202; lLanzetta. Turnshek fc Sandoval 
1993) do not show such a strong Lyman edge absorption. 

Alternatively, our Balmer edge absorption feature could 
be imprinted in the scattering medium itself, i.e. both scat- 
tering and absorption are occurring in the same region. This 
is actually the favored model for the circumstellar disks of Be 
stars. In these systems, the radiation from the stellar photo- 
sphere is electron-scattered by the circumstellar disk. Due to 
the absorption opacity in the disk medium competing with 
electron scattering, absorption signatures are imprinted in 
the scattered light spectrum. Thus a large Balmer ed ge is 
seen in the polarized flux (e.g. iBiorkman et alJ Il99ll note 
that their Fig.l shows Fx and P, not directly P x Fa). 

The same thing might be happening to our quasars. 
The edge feature could be imprinted by the circumnuclear 
scattering region — this scattering region might possibly 
be the BLR clouds. However, our argument is based on the 
empirical fact that there is essentially no feature seen in the 
broad-line wavelengths, either in absorption or in emission. 
This fact (1) suggests that the spatial scale of the scattering 
region is smaller than that of the BLR as we have argued, (2) 
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seems to be against any simple scattering model involving 
in situ absorption. 

Also, the fact that the edge feature looks quite broad- 
ened might require quite specific conditions (e.g. high veloc- 
ity dispersion) for the medium producing the feature. If the 
feature is imprinted in the scattering region, these conditions 
might put the scattering region as deep in the potential well 
as a part of the continuum emitter: this is equivalent to one 
of our simple pictures, i.e. the polarization and the feature 
in the polarized flux originate from the BBB emitter and 
thus are intrinsic to the BBB emission. 

It is possible that the outer part of the disk itself scat- 
ters the radiation from the inner regions. This would give 
the right polarization direction parallel to the disk symme- 
try axis (which is supposed to be along the radio jet axis). 
The Balmer edge feature could be intrinsic to the radiation 
from the inner regions or could be imprinted by the outer 
part of the disk. In either case, the medium producing the 
feature must satisfy the constraints we discussed above: that 
the broad lines are essentially unseen and the edge is broad- 
ened. 



4.2 UV spectral features 

We have so far concentrated on the spectral feature at the 
Balmer edge region. However, there are apparently more 
spectral features in the polarized flux at shorter wavelengths 
in at least some of the objects shown in Fig l31l a local peak 
at ~ 3050 A, and/or absorption shortward of ~ 3050 A, and 
a possible absorption at ~ 2600A. These look clearer in the 
overplotted and composite spectra shown in Fig l34l and 1351 
To produce the composite, we first normalized the five polar- 
ized flux spectra by the mean flux between 2800A and 5000A 
in the rest frame of each quasar excluding the wavelengths 
affected by telluric absorptions, binned the spectra with a 
20A width (shown in the upper panel), and then took the 
average of these five with equal weighting (the lower panel) . 

The identification and interpretation of these possible 
spectral features is not still clear, but they might be re- 
lated to Fell absorption and possibly the Bowen resonance- 
fluorescence lines which are expected often to be highly po- 
larized. One of the authors (RA) would offer US$100,000 
to anyone who successfully identifies them. To be decisive 
on these possible features, we need more data with a higher 
S/N at this wavelength region on these and other quasars. 

4.3 Radio properties 

In Table [5] we have summarized the radio properties of our 
quasar sample along with the polarimetric properties. The 
quasars observed are all radio-loud except Q0205+024. For 
these radio-loud quasars, radio maps with resolved struc- 
tures are available in the literature except for Q1912-550. 
All of them show edge-brightened FR-II morphology. 

It is known that there is a strong statistical ten- 
dency for the optical polarization E-vector direction to 
lie parallel to the radio structural axis in quasars and 
Seyfert Is JStockman et alJll97gt jMoore fc Stockmanlll984t 
iRusk fc Seaauisdll985tlAntonuccll983l.l2002l) . though there 
are some cases where the optical polar ization direction lies 
rather perpendicular to the radio axis (M oore fc Stockman! 
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Figure 36. The histogram for the difference of the radio axis PA 
and optical polarization PA. 

I1984 lAntonucdl2002t ISmith et al.ll2002l) . Some quasars in 
our sample also show a perpendicular relationship (2 of them 
are included in Moore and Stockman sample), while most 
of the others show parallel polarization. In Fig l36l the his- 
togram of the difference of these two PAs (APA) is shown 
for our sample. It is not certain whe ther this distribution 
indic ates bimodality (e.g. a dip test jHartigan fc Hartiganl 
1985] does not reject unimodality with a significant confi- 
dence level). The bimodali ty is also uncertain i n the large 
samples such as those of iMoore fc Stockman! lll984l) and 
IRusk fc Seaauisd <ll985l) (again, a dip test does not give a 
statistically significant result). In Tabled we have classified 
the quasars into three groups simply as 'parallel' (APA < 
30°), 'intermediate' (30° < APA < 60°), and 'perpendicu- 
lar' (60° < APA). Some properties seem to have a possible 
APA dependence as we discuss below. 

4-3.1 rp, line polarization and Balmer edge 

Most of the quasars in the parallel group show a decrease 
in polarization shortward of 4000A, meaning an (at least 
partially) unpolarized small blue bump. We show the polar- 
ization ratio rp in Table [5] (taken from Table |SJ . Those in 
the parallel group also show essentially unpolarized emission 
lines: many of them have a tight upper limit on the broad 
line EW in the polarized flux compared to that in the total 
flux, as we show in Table |U] the ratio of the former to the 
latter (for H/3, or Mgll if H/3 not available; see Tables|7|and 
[SJ. It also seems that the quasars in the parallel group tend 
to have the Balmer edge absorption feature in the polarized 
flux: all the five quasars shown in Fig l31l are in this parallel 
group, and the feature is possibly seen in another parallel 
quasar Q0003+158. 

On the other hand, the quasars in the perpendicular 
group do not seem to show the P decrease shortward of 
4000A (r P ~ 1 or > 1), which indicates that they tend 
to have a polarized small blue bump. While the constraint 
on rp is still not tight in two of the three quasars, the 
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Table 9. Radio properties along with polarization properties. Luminosities are calculated assuming Ho = 70km/s/Mpc, f2 m = 0.3 and 
Ca = 0.7. The values for uL v are from our spectroscopic measurements, so they are only approximate. L§q is the radio luminosity at 
rest-frame 5GHz, calculated assuming /„ oc u~ 0,5 . PA op t and PA cor are the PA of t he optical polarization at 40 00-473lA before and after 
the IS P c orrection, respectiv ely (from Table |tj). Referenc e fo r the radio map is: 1 . . iHutchings et all ll99Sj). 2jReWjKronb£r£^&Perle\J 
ll99Stl . 3.. lPrice et al] il993ft . 4.. lKellermann et alj il994D . 5.. lRudnick. Sitko. fc Steinl il984h . 6.. lAntonuccil il985l) . 7.. I Fanti et al .1 Jl977fl 
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other one (Q2115-305) shows P increase instead of decrease. 
The broad lines also look polarized in these objects: this 
is quite certain in Q2115-305; marginal in Q0350-073; in 
Q2251+113, the line polarization is puzzling in that the 
PA rotation seen at the lines is not present in the small 
blue bump wavelengths. We have plotted rp against APA in 
Fig l37l which illustrate the correlation between them (95.9% 
statistical significance in the Spearman's rank correlation 
coefficient, and 96.2% in the Kendall's). 

Note also that rp depends on the strength of small blue 
bump in the total flux. In fact, the small blue bump seems 
to be stronger in the parallel ones and weaker in the perpen- 
dicular ones. This is shown in Fig l38l As a rough strength 
measure for the small blue bump, we simply take the ra- 
tio of the total flux at 3050A to that at 4200A (so of course 
this also depends on the color of the underlying continuum) , 
and convert it to the corresponding spectral index a where 
/„ oc v a . This index is plotted against APA, which illus- 
trates the correlation (99.8% statistical significance in the 
Spearman's rank correlation coefficient, and 99.9% in the 
Kendall's). Thus, r p is probably measuring the effects from 
both the polarization and the strength of the small blue 
bump. 

For the intermediate group, we note that Q0405-123, 
as we described in the previous section, shows a polariza- 
tion PA rotation shortward of rest 4000A, which is similar 
to 4C09.72 (and 3C95 or 3C323.1) and may be due to ISP 
contamination. The PA of the intrinsic polarization could be 
at ~160°, which would put this object close to the 'parallel' 
group (APA ~ 33°). Note also that Q1004±130 may have an 
ISP contamination, and the intrinsic polarization PA could 
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Figure 37. The ratio of the polarization at 2891-3600A to that 
at 4000-4731A is plotted against APA. 



be at ~45° , which would put it into a 'perpendicular' group 
(APA ~ 76°). Even before this correction, this object looks 
similar to Q2115-305 in terms of the rising P toward shorter 
wavelengths. These possible re-evaluations of APA look con- 
sistent with the properties described above for the parallel 
and perpendicular groups, though we did not detect H/3 line 
polarization in Q1004+130. 
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Figure 38. The spectral index a (/„ oc u a ) corresponding to the 
flux ratio at 3050A and 4200A is plotted against APA. 



4-3.2 Rv and viewing angle 

For the ones with resolved radio structure, the ratio R 
of the core to lobe luminosity can provide a rough in- 
dication for the viewi ng angle, the angle o f our line of 
sight to the jet axis ijWills fc Browne! Il98fift . Or, a bet- 
ter quantity is the ratio of the core lu minosity to opti- 
cal luminosity J Wills fc Brothertodll995h . denoted as Rv. 
We have calculated this latter ratio Rv for each quasar in 
our sample with ava ilable core flux measurements, following 
IWills fc Brothertonl dl995T) . The values are shown in Tabled 
Note that the value for Q2115-305 is regarded as an upper 
limit since its available map is coarse and no flat-spectrum 
core seems to be well isolated. The values for the two per- 
pendicular ones and Q1004+130 are smaller than the others, 
suggesting that the perpendicular ones have a larger view- 
ing angle than the parallel ones (this tendency is also seen 
in the values of R). 



4-3.3 Possible interpretation for these and other AGN 

The overall tendencies for the small blue bump polarization, 
broad line polarization, and viewing angles could be under- 
stood if we suppose that (1) electron scattering is occuring 
in a flat or oblate compact (~100-1000r 9 ?) optically thin re- 
gion interior to the BLR, and (2) electron or dust scattering 
is occuring in an extended (> pc) optically thin polar region 
exterior to the BLR, as deduced from the unified model. 

When the viewing angle is small, the scattered light 
from the ~100-1000r 9 flat compact region dominates, result- 
ing in a small parallel polarization of the continuum only. 
The objects have Type 1 spectra. When the viewing angle is 
somewhat larger, and the the inner disk-like flat scattering 
region becomes occulted or otherwise ineffective, the >pc- 
scale region dominates: this would cause a small perpendic- 
ular polarization when the BBB and BLR aren't occulted or 
just attenuated leading to a perpendicular-polarized Type 1, 
and a Type 2 when these components do become occulted. 
In fact the continuum of the perpendicular ones in our sam- 



ple might be redder than the parallel ones, based on the 
flux ratio at 3050A and 4200A which would depend both on 
the color of the continuum under the small blue bump and 
the strength of the small blue bump (the continuum color 
longward of the small blue bump wavelengths could not be 
obtained reliably from our data for high-redshift ones in par- 
ticular). 

For the polarization of Seyfert 1 galaxies, Smith et al. 
(2004) have already suggested that two different scattering 
regions exist, one in an inner equatorial region and the other 
in an outer polar region. The difference from our case is that 
those Seyfert 1 galaxies with the parallel property have po- 
larized broad lines, and the polarization PA is different from 
that of the continuum with the p olarization magnitu de lower 
than that of the continuum fe.g. lSmith et al] |2002j^ This is 
also t he case in some broad-line radio galaxies (ICoherTet alJ 
1999). In these cases, the inferred inner equatorial scattering 
region would be comparable in size to the BLR. 

The distinction between the ones with essentially unpo- 
larized lines (i.e. most of the parallel quasars in our sample) 
and the ones with partially polarized lines is not yet clear, 
but it might be related to optical luminosity. Many of the 
objects with partially polarized lines are of a rather low op- 
tical luminosity [yL v at 4000A < 10 45 erg/s), while the ones 
with unpolarized lines in our sample are of a rather high op- 
tical luminosity (yL v at 4000A > 10 45 erg/s; see Table 0. 
However, an object like 3C382, which has partially polar- 
ized broad lines, has a rather high optical luminosity {vL v 
at 4000A ~ 2 x 10 45 erg/s). We need a larger sample for the 
distinction to be clarified. 



5 CONCLUSIONS 

In our previous paper, we have reported our detection 
of a Balmer edge feature in one quasar, Ton202, which 
was observed in May 2002. We have since implemented 
VLT/FORS1 and Keck spectropolarimetry in Sep 2002 and 
May 2003, respectively, on 14 more quasars with a known 
~1% optical polarization, and also re-observed Ton 202 in 
May 2003. Some of these observations are affected by inter- 
stellar polarization, but the measurements have been cor- 
rected for this effect reasonably well. Among the newly ob- 
served 14 quasars, the high or relatively high S/N polarized 
flux of at least four quasars (4C09.72, 3C95, 3C323.1, B2 
1208+32) shows essentially no emission lines. We have found 
that these four quasars show a Balmer edge feature in the 
polarized flux, which is quite similar to that seen in Ton202 
observed in May 2002. 

Based on the unpolarized broad emission lines, the po- 
larized flux is thought to originate interior to the BLR, and 
the Balmer edge feature is most simply interpreted as an 
intrinsic spectral feature of the Big Blue Bump emission. 
In this case, the edge absorption feature indicates that the 
BBB emission is thermal and optically thick. However, we 
have also discussed other alternative interpretations. 

The Balmer edge feature identified in the polarized flux 
of these five quasars is entirely a new finding. The feature 
is apparently widespread, based on its similarity in these 
well observed objects, and the possible similar behavior of 
the polarized flux at around 4000A in two other objects 
(Q0003+158, Q1912-550). There also seem to be more spec- 
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tral features in the shorter rest UV wavelengths of the po- 
larized flux, which are common to multiple objects. 

In the re-observation of Ton202 in May 2003, we did 
not find a similar dramatic feature. This is apparently due 
to real polarization variability. This is not unexpected since 
the polarization is thought to originate from a small scale, 
less than the size scale of broad line region. The polariza- 
tion variability of normal quasars (aside from synchrotron 
dominated objects) has been reported in the literature. This 
object, as well as other objects should be followed up to con- 
firm the variability. 

Among the 15 quasars in our sample, 13 are radio-loud 
quasars with an available radio map. While many of them 
show the optical polarization roughly parallel to the radio 
structural axis, some of them show the polarization rather 
perpendicular to the radio axis (however, the bimodality 
is not certain). Although the sample is too small to draw 
strong conclusions, the parallel ones tend to show an unpo- 
larized and stronger small blue bump (thus shows a polar- 
ization drop shortward of 4000A), unpolarized broad lines, a 
Balmer edge absorption feature in the polarized flux, and a 
smaller viewing angle, while the perpendicular ones tend to 
show a polarized and/or weaker small blue bump, polarized 
broad lines, and a larger viewing angle. 

Based on our spectropolarimetric data and other data 
in the literature, the parallel type 1 objects seem to have 
at least partially unpolarized broad lines. The ones with 
(almost completely) unpolarized lines might have a higher 
optical luminosity than the ones with partially unpolarized 
lines. However, a much larger sample is required to clarify 
this distinction. 
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